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Abstract 

We study information geometry of the Dicke model, in the thermodynamic 
limit. The scalar curvature R of the Riemannian metric tensor induced on 
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1 Introduction 



The physics of phase transitions has been the focus of intense research over the 
last century. Information theoretic analyses of the same are relatively recent, 
and have mostly found applications in relation to second order classical phase 
transitions (CPTs) pQ. The essential idea herein is geometric in nature [2] and 
rests on the definition of a Riemannian metric on the equilibrium thermodynamic 
state space of a system. Ruppeiner [3] has conjectured that the scalar curvature 
R arising out of such a metric is related to the correlation volume of the system, 
i.e \R\ ~ £ d , where d is the system dimension. This association follows from 
thermal fluctuation theory, with the understanding being that the greater is the 
"distance" between two equilibrium thermodynamic states, the lesser is the prob- 
ability that these are related by a fluctuation (for a nice, pedagogical explanation, 
see [4]). The relation between the scalar curvature and the correlation volume has 
been tested for a variety of classical models [3j and has also been applied to the 
theory of first order phase transitions in [5], and show excellent agreement with 
experimental results for liquid-gas systems. (For related recent results in liquid- 
liquid phase transitions and for Lennard- Jones fluids, see [6], and [7]). It is by 
now established that for classical systems, R diverges on the spinodal curve, and 
its equality in co-existing phases indicates a first order phase transition. These 
two aspects can be combined to map the phase diagram of the system, purely 
from a geometric perspective. 

Information theoretic methods can be applied to quantum mechanical sys- 
tems as well. The pioneering work in this direction was carried out by Provost 
and Vallee [Sj, who defined a Riemannian metric tensor from the Hilbert space 
structure of quantum states of a system. The application of this method to the 
study of zero temperature second order quantum phase transitions (QPTs) is 
more recent, and appeared (following the results of |9J) in the seminal work of 
Zanardi et. al. [10]. In particular, these authors showed that the scalar curvature 
arising out of the metric on the manifold of coupling constants^ in the transverse 
XY spin chain is an indicator of QPTs in the system. 

Let us remind the reader that information geometry of QPTs has very different 
qualitative features compared to that of CPTs. This is not unexpected, given 
that the physical mechanisms involved in the two situations are different. For 
classical systems, the geometric parameters are thermodynamic variables and 
their Legendre transforms (hence information geometry in the context of CPTs 
is also referred to as thermodynamic geometry), whereas for quantum systems, 
the geometry is characterized by parameters that appear in the Hamiltonian, i.e 
the coupling constants of the theory. The relation between the scalar curvature 
and the correlation length (\R\ ~ £ d ) that holds for CPTs, does not, in general, 
hold for QPTs, as we will show. 

The work of [10] puts information theoretic studies of QPTs into perspective, 
along with the more popular methods involving fidelity, defined as the modulus of 
the overlap between two infinitesimally separated ground state wave functions [9], 
and fidelity susceptibility (see, e.g. pj]. For an excellent review, see [12] )■ How- 
ever, although the literature on the fidelity approach to the study of QPTs is 

lr TIiis will be referred to as the parameter manifold, or the information geometric manifold. 
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vast, less attention has been paid to the geometric structure of the parameter 
manifold, and in particular, its scalar curvature R. While for classical spin mod- 
els and fluid systems, R always diverges at a second order phase transition, the 
same is not true for QPTs. For example, it was shown in pi)] that for the trans- 
verse XY spin chain, R develops a discontinuity at the Ising transition lines and 
does not diverge there. While an argument arising out of first order perturba- 
tion theory shows that the components of the metric tensor on the parameter 
manifold generally diverge at a QPT [10] (this may not be true in some specific 
examples, see [H]), the same is thus not true for the scalar curvature. The latter 
is however a useful quantity to study, since it is invariant under coordinate trans- 
formations (of the parameter space), as opposed to the metric components (that 
transform as rank two tensors under such transformations). For two dimensional 
parameter manifolds which is the case we will be concerned with in this paper, 
standard textbook results indicate that the Riemann curvature tensor |14J has 
only one independent component, and hence leads to a unique scalar curvature. 
Whereas divergences in the components of the metric tensor might be due to 
specific coordinate choices, and can be removed by (possibly singular) coordinate 
transformations, those of the scalar curvature correspond to ones that cannot 
be removed by such transformations, and therefore characterize the global (i.e, 
coordinate independent) properties of the manifold. 

It is important and interesting to study the behavior of the parameter mani- 
folds, for models exhibiting QPTs, via its scalar curvature, as this is a coordinate 
free approach and is further expected to give insights into the difference in the 
geometrical nature of QPTs vis a vis CPTs. To the best of our knowledge, this 
has not been done before. The main purpose of this paper is to study the be- 
havior of R for the Dicke model of quantum optics, in the thermodynamic limit 
(TDL), for finite values of de-tuning. We show that in this model, although the 
components of the metric tensor diverge at QPT, R is regular at, and continuous 
across, the phase boundary, and hence the parameter manifold is smooth, even at 
the location of the phase transition. Our results indicate that the global behavior 
of the parameter manifold varies widely between models, and does not allow for 
an universal characterization of QPTs, a situation that is markedly different from 
the classical case. 

This paper is organized as follows. In the next section, we briefly review 
information geometry and the results of pU] for the transverse field XY model. 
This section is mostly meant to set the notations and conventions for what is 
presented in the rest of the paper, but towards the end, we indicate some novel 
features of R. In section 3, we analyze the Dicke model of quantum optics, and 
calculate the associated scalar curvature, and the fidelity susceptibility. This is 
done both with and without the rotating wave approximation. Section 4 ends 
with our conclusions and directions for future study, where we also indicate some 
preliminary results on similar analyses for the Lipkin-Meshkov-Glick model and 
the Kitaev honeycomb model. 
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2 Information Geometry and the Transverse XY 
Model 



In this section, we briefly review the basics of information geometry, and the 
results of [10] in the context of the transverse field XY model. The essential idea 
here [8] is to use two infmitesimally separated quantum states, and construct 
(upto second order) the quantity 



\ip (q + dq) — tjj (q) | 2 = (cV0|c^) dq 1 dq j = cxijdq l dq j 



(1) 



where q % (collectively denoted as q in the l.h.s of eq.((T])) denotes the parameters 
on which the wave function depends on, and di is a derivative with respect to 
q % . The a^s are not gauge invariant, and in [8], it was shown that a meaningful 
second rank tensor can be defined as 



9ij = a ij - A = -ity (?) \9iip (<?)) 



(2) 



where the quantities ^ specify a metric on the parameter manifold, induced 
from the natural structure of the Hilbert space of quantum states. This metric 
can be shown to be positive definite (i.e the diagonal elements and the principal 
minors are positive definite), and hence one obtains a Riemannian structure via 
a distance, i.e a line element ds 2 = gijdq % dqi which can be used to study the 
system under a variation of the parameters q % . Note that the physical meaning 
of this metric (and the line element) lies in the fact that the larger the distance, 
the higher is the statistical distinguishability of the two states. For the purpose 
of studying quantum phase transitions, it is natural to focus on the ground state 
overlap and the corresponding induced metric [9], [10] . For the case of a two- 
parameter model, i.e q % = (q l , q 2 ), the scalar curvature for the manifold with this 
metric can be computed via the standard formula R 



—j=(A+B), where 



A 
B 



_d_ 

dq 1 
_d_ 

dq 2 



9i2 dg 



ii 



1 dg: 



22 



_0n Va d( l 2 V9 dq 1 

2 dg 12 1 dgn 

y/g dq 1 y/g dq 2 g lly /g dq 1 



9i2 dg 



ii 



(3) 



Here, g^-, i,j = 1,2 are a function of the the parameters (g 1 ,? 2 ), and g is the 
determinant of the metric of eq.([2]). 

One of the few models that yield analytic results for the scalar curvature is 
the transverse XY model [TO], which we briefly review. The Hamiltonian for the 
model with (2iV + 1) spins is 



H 



XY 



N 

E 



1 + 7 



a 3 a j+1 + 



1 



-N 



1 y y h z 



(4) 



where the a 1 



x, y, z are Pauli matrices, 7 is an anisotropy parameter and h 



is the magnetic field. Diagonalizing this Hamiltonian with the Jordan- Wigner, 
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Fourier and Bogoliubov transformations, it can be checked [ID] that with the co- 
ordinates (q 1 , q 2 ) = (h, 7), the components of the metric induced on the parameter 

manifold is given by gfy = \ J2k>o ft' wnere = {(h - Cos(x k )) 2 + 7 2 Sin 2 (x fc )) 2 
and 

hh = 7 2 Sin 2 (x fc ); # 77 = Sin 2 (a; fc ) [h - Cos(x fc )] 2 ; g hl = 7Sin 2 (a; fc ) [Cos(x fc ) - h] 

(5) 

where x k = L = 2N + 1. In the TDL, the sum over k can be replaced by an 
integral, which can be evaluated in the complex plane [10], and it can be shown 
that metric components have four poles (only g 77 has an additional pole at the 
origin), located on the complex plane at values of the coordinates 

± _ h± y/fr 2 + 7 2 -l , ± _ h± V /h*+ 1 *-l 

Al " T^f ; As " 1 + 7 (6) 

Evaluating the integrals then imply that some of the metric components are 
divergent at 7 = and at h = ±1 (see below), and the scalar curvature is 

32\ 1 

TJ r~Y;:.. 

32\ 1 



R = -[ — ) r+T^' \h\<l, J>0 



\h\ > 1, 7> (7) 



L J 1 - A^A 2 

with similar expressions for 7 < 0. This is equivalent to the result derived in 
(as can be checked by substituting the values of \f 2 in eq.([7|)), but the usefulness 
of eq.([I]) is that, using the results of [IE], it reveals that for the region \h\ < 1, 
\R\ ex (l — e -1 /^) , where ^ is the correlation length inside the circle /i 2 + 7 2 = 1. 
No such relation exists for the region \h\ > 1. Hence, on the line 7 = (with 
\h\ < 1), the scalar curvature diverges, as does the correlation length, but this 
is not the case for the Ising transition lines. This indicates that, as alluded to 
in the introduction, there is no simple relation between the scalar curvature and 
the correlation length at a second order quantum phase transition, as opposed to 
classical cases. 

We also emphasize that in this example, the scalar curvature R does not show 
any special behavior near the Ising critical lines. It does show a discontinuity 
along these lines [ID], but no geometrical meaning can be attributed to it, since 
the information metrics are different in the regions separated by these lines. The 
only genuine (i.e coordinate independent) curvature singularity is on the line 
7 = for \h\ < 1, where R diverges, and is related here to the correlation length, 
as we have mentioned. The singularities of the metric near the Ising transition 
lines are thus simply artifacts of coordinate choices. For example, in the region 
\h\ < 1, the metric can be written in two equivalent forms, 



^ = -(4. + ^)^^ + ^) (8) 
16 7 Vl-^ 2 (1 + 7)/ !67V (1 + 7) 
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by a redefinition of coordinates h = Sin -1 (h). Whereas the first metric of eq.® 
has a divergence in the limit \h\ — > 1 in its first component, the second metric 
does not have any divergence (apart from that at 7 = which is a genuine 
singularity and cannot be removed by a coordinate re-definition). The coordinate 
transformation indicated above is undefined at h = ±1, but note that in the 
region \h\ < 1, it gives rise to the same scalar curvature as does the first metric of 
eq.©, and thus describes globally the same manifold as the first metric. For the 
region \h\ > 1, the metric can be calculated in a similar fashion, and it turns out 
that R does not show any divergence in this region as well, except at the multi 
critical point (MCP), 7 = 0, \h\ — > 1. The information geometric manifold, in 
this region, is thus smooth (excepting for the MCP), and the divergences of the 
metric tensor can again be removed by appropriate coordinate transformations. 



3 The Dicke model with a single Bosonic mode 

We now use the formalism described in the previous section to consider the in- 
formation geometry of the Dicke model (see, e.g [17], [18], [IS]), which is of great 
relevance in quantum optics. This model describes a bosonic mode interacting 
with a system of N two-state atoms, via a dipole interaction, and the Hamiltonian 
is given by (we will set h = 1 in what follows) [20J 

H = oj J z + ujata + (a f + a) ( J + + J J) (9) 

Here, a and a) represent the bosonic mode annihilation and creation operators 
respectively (of frequency u). The J's are the collective angular momentum 
operators obtained by summing over the atomic spin operators, and A is the atom- 
field coupling strength, which drives QPTs in the model, and will be referred to 
as the driving parameter of the model in sequel. Uq is the difference between the 
energies of the two states of the atom. We also set j = —, where N originates from 
the volume factor of the cavity in the original Dicke Hamiltonian [17] . In [20] . 
it was shown that this choice of j describes the iV-atom two-level system by a 
single (N + l)-level system. 

In order to understand the information geometry of the Dicke model, we 
treat both u> and A as geometric parameters. Effectively, we consider a family of 
Hamiltonians of the form in eq.(j9]), which smoothly depends on u and A. From 
a geometric perspective, it is natural that all the parameters of the model are 
treated at par, and we are simply looking at a class of models with non-zero 
values of de-tuning. 

We will first consider the Dicke model in the TDL, using the rotating wave 
approximation (RWA). The Hamiltonian is given by 

Hrwa = UoJz + (jja)a H -=(a) + aJ + ) (10) 

To diagonalize this Hamiltonian in the normal phased we follow the method of [20] . 
First, we record the standard Holstein-Primakoff representation of the angular 

2 The presence of extra terms in the Hamiltonian for the super-radiant phase in the RWA 
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momentum operators 

J + = 6V 2 J - ^b; J- = ~ tfb b; J z = tfb - j 
and in the thermodynamic limit N — > oo, we obtain, using these, 



Using the position and momentum operators 

1 



3^0 



x 



y 



a f + a); p x = l \j^{ a> - (, ' ! 



the Hamiltonian of eq. (|T2|) is given by 
1 



RWA 



w x + p x + co y +P y + 2\^tuuj xy + 



U! — U 



(12) 



(13) 



-J wo (14) 



Now, we rotate the coordinate and momentum axes by 



" R ' 




Zi 


o" 
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P _ 







z 2 _ 







(15) 



where is the 2x2 null matrix, and R = (x,y) T , P = (p x ,P y ) T , q = (91,92) 
p = (px,p2) ■ Zi and X2 are 2x2 matrices given by 



T _ ( C0S7 Sin 7 \ q- _ ( Cos 7 ' ShryM 
Ll ~ ^-Sin 7 Cos 7 J ' X2 - \ -SinV Cosy J 



(16) 



Diagonalizing the Hamiltonian of eq. (H4j) then sets Tan2 7 = an d l' ~ \ ■ 



and we have two uncoupled oscillators : 



"0 



RWA 



2 2 1 2 1 1 

£-5i +Pi 1 



A 



+ e 2 + q 2 2 +pl I 1 



A 



•- (w + w )-ju;o 
(17) 



where we have defined the collective excitation modes as e± (denoted as the 
atomic and photonic branches in [20] due to their values at A = 0) : 

e± 2 = \ [ul + u 2 ± sJ{uI-^Y + A^uju^ (18) 



makes it difficult to diagonalize it using position and momentum representation, and we will 
confine ourselves to the normal phase in this analysis. 
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The Hamiltonian of eq. (fT7|) is then re-quantized with two new bosonic modes, 



(c\ - Cl ) 

(4 - Pa) (19) 

- (u Q +u)- ju (20) 

where m± = 1 ± X/^/uuo, and e± = e±/ ■Jm±. The ground state wave function 
for the model can now be obtained as 

M*, y) = f ^'^jr ^ 4 " -) V4 e- (1/2)(M > (21) 

and the overlap between two infinitesimally separated ground states is 

(W ;> = 2 ' D ^ DeU 'i! (22) 
[Det (A + A')] 1 

where A = U^MU, with 

M = diag( e „ v ^7, e+V ^T); U = ( ^ (23) 

and A' is obtained by infinitesimally shifting A — >• A' and u — > u' . 

Given the overlap of the ground state wave function of eq. fl22|) . one can cal- 
culate the components of the metric induced on the parameter manifold by using 
eq.([2]). The reader will notice that obtaining analytical results for the same is 
difficult, and we will present numerical results. First, note that in this case, there 
is a quantum phase transition at A c = y/cou [20] where the excitation energy of 
the photonic branch, e_ = 0. Setting uj q = 1 for convenience, we find that all 
components of the metric tensor diverge at A c = y/uj, as expected. We point out 
that the diagonal components of the metric tensor diverge to oo, whereas the off- 
diagonal components diverge to — oo, and the determinant of the metric tensor 
diverges to positive infinity. However, we find that there is no divergence of the 
scalar curvature at the phase boundary. This indicates that the underlying man- 
ifold in the parameter space does not exhibit any special behavior here, and we 
conclude that the divergences of the components of the metric tensor are artifacts 
of the coordinate choice, as was the case of the transverse XY model of the pre- 
vious section (along the Ising transition line). We find that the scalar curvature 
R — > 4 at the phase boundary for all values of u. We also calculate the fidelity 



1 



2e' 



with its final form being 



c\ + cx) , pi= i\l — 




4 + c 2 I , p 2 = i 




Hrwa = m_e'_ ^4 c i + +m + e' + (c\c 2 + 
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Figure 1: (Color online) The scalar cur- 
vature R of the Dicke model in the RWA 
with Uq — 1, plotted function of 
the driving parameter A in the normal 
phase, for oj = 0.6 (solid green), uj = 0.8 
(dashed blue) and uj — 1.2 (dot-dashed 
orange) . 



Figure 2: (Color online) Log of fidelity 
susceptibility cLS cL function of A in the 
RWA for the Dicke model in the normal 
phase with uj — 1, for the paths uj = r/X, 
with r] = 1 (solid green), 2 (dotted blue), 
3 (dot-dashed orange) and 4 (dashed 
black) . 



susceptibility along the paths parametrized by uj = 77 A for 77 = 1, • • - 4 using the 
methods of [12] . We find that, being a sum of the metric tensor components, 
this diverges at the phase boundary as expected, although this divergence does 
not have any bearing on the global structure of the parameter manifold. Our 
results are shown in figs. ([I]) and (|2]), where we have shown numerical plots for 
the scalar curvature R and the log of the fidelity susceptibility, as a function of 
the driving parameter A, for different values of uj. We note that the peaks in the 
scalar curvature in fig. (JTJ) are possible mathematical artifacts and do not carry 
any physical significance. They are absent when we consider the theory without 
the RWA, to which we presently turn to. 

For the Dicke model in the TDL, without the RWA, the results are qualita- 
tively similar. Here, the overlap between two infinitesimally separated ground 
states in the normal phase El is given by [9], [20] . 



("00 (w, A ) I -00 A')) = 2 



DetA Bet A' 



Det [A + A' 



(24) 



where the matrix A is 



.4 



l/ / «-/3Cos27 /3Sin27 \ 
2 V /3Sin27 a + /3Cos27 J 



(25) 



3 All the quantities are denoted by a tilde here compared to the analysis with the RWA, to 
avoid confusion. 
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Figure 3: (Color online) The scalar cur- 
vature R of the Dicke model without the 
RWA, with ujo — 1, plotted as a function 
of A, for uj = 0.8 (solid green), Cj = 0.999 
(dotted blue) and Cj = 1.2 (dot-dashed 
orange) . 



Figure 4: (Color online) Log of the Fi- 
delity Susceptibility for the Dicke model 
without the RWA, for ojq — 1, plotted as 
a function of A, for the paths Cj = A (solid 
green), Cj = 1.5A (dotted blue), Cj = 2A 
(dot-dashed orange). 



with a,/3 = e+ ± e_, 7 = |Tan _1 f 4 ^p° j- The matrix A' is obtained by 

replacing a) — > Cj' and A — Y X' in the matrix A, and similarly for e±. Here, 
e± are the energies of the two uncoupled oscillator modes (atomic and photonic 
branches), given by 

&± = \ (&o + ^ ± \/(w 2 -w 2 ) 2 + 16A 2 u;u;o) (26) 



At A c = |-y/SSo, e_ = 0, and the system undergoes a QPT. 

For the overlap function of eq.( l2~4l) . one can calculate the metric components 
and the scalar curvature just as before, by treating Cj and A as geometric param- 
eters. Again we find that the diagonal (off-diagonal) components of the metric 
tensor diverge to 00 (—00) and the determinant of the metric diverges to 00. A 
similar analysis can be carried out in the super-radiant phase using its ground 
state wave function derived in [20], and the scalar curvature R can again be 
obtained as a function of A and Cj (we set Cjq = 1). Combining this with the 
analysis for the normal phase, we obtain R in the entire region of the parameter 
space. In fig.()3]), we plot R, as a function of the driving parameter A for various 
values of Cj. In the normal phase, the scalar curvature does not show any peak 
(present with the RWA) and its value remains ~ 4 throughout, consistent with 
what we obtained with the RWA. In the super-radiant phase, R — > 4 at the phase 
boundary, and shows a maximum, slightly away from the boundary, but we are 
not aware of a physical reason for the same. We also find that in the limit of 
large A, R — > 0, which is consistent with the fact that in this limit, the modes 
e + and e_ decouple [20], and the ground state wave-function is independent of 
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A. As a result of this, the parameter manifold becomes one-dimensional in this 
limit, and hence has zero curvature. R being a coordinate independent quantity, 
our conclusion here is that the information geometric manifold is smooth at the 
location of the QPT. The scalar curvature is also seen to be continuous across 
the transition, from fig.(|3]), in contrast to the result for the transverse XY model. 
Finally, fig.flU) shows the results for the log of the fidelity susceptibility, which 
expectedly diverges at the phase transition. 



4 Discussions and Conclusions 

In this concluding section, for the sake of completeness, let us first mention some 
cases that we believe, deserve further attention and analysis. Our first example 
is that of the Lipkin-Meshkov-Glick model for N spins (see e.g [21] for a finite 
temperature analysis, and references therein), given by the Hamiltonian 

(*<i) » 

where the ex's are the Pauli matrices, and 7 and h are the anisotropy parameter 
and the magnetic field respectively. After diagonalizing the Hamiltonian using 
Holstein-Primakoff and Bogoliubov transformations (we omit the details of the 
standard calculations here), it can be shown that the induced metric on the 
parameter manifold has the components 

(7 -I) 2 1 7-1 (l?R s 

9hh 32(/i-l) 2 ( 7 -/i) 2 ' 911 32(>y - h) 2 ' 9hnf 32(/i-l)( 7 -/i) 2 1 ' 

It can be seen that the determinant of the metric is zero, signaling that informa- 
tion geometry is ill defined in this model. A similar conclusion has been reached 
by different methods in [22], and this needs further investigation. In the context 
of the Kitaev honeycomb model, we find that the parameter manifold is smooth, 
i.e does not show any divergence. The scalar curvature here shows the same 
qualitative behavior as the fidelity susceptibility [23], i.e R show a series of peaks 
in the gapless region, exactly coinciding with the peaks of the fidelity suscep- 
tibility [23J. Information geometry of this model also needs to be understood 
further. 

To summarize, in this paper, we have undertaken a critical analysis of infor- 
mation geometry in quantum phase transitions, in the context of the Dicke model 
of quantum optics in the thermodynamic limit. By analyzing the scalar curvature 
of the parameter manifold, which relates to coordinate independent properties of 
the same, we have shown that this manifold is smooth for the Dicke model, and 
does not show any singular behavior at phase transitions, although the compo- 
nents of the metric tensor diverge there. We have studied the geometry for the 
Dicke model, both with and without the rotating wave approximation, and found 
consistent results for the scalar curvature of the parameter manifold. We have 
also shown that the scalar curvature is continuous across the phase boundary. 
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Our results indicate that purely from a geometric perspective, the divergences 
in the metric tensor in this model can be removed by appropriate coordinate 
transformations. 

Our results also point to the fact that information geometry does not seem to 
provide a universal characterization of second order QPTs via the scalar curvature 
of the parameter manifold. Apart from the Dicke model, where R is regular, we 
have demonstrated this explicitly in the context of the XY spin chain, where it 
was shown that R, in general, does not have a simple relation with the correlation 
length at a QPT. These statements, along with our discussion in the beginning 
of this section on the LMG and the Kitaev honeycomb models, indicate that for 
QPTs, the geometry of the parameter manifold needs to be studied on a case by 
case basis. 

This should be compared to the case of classical phase transitions, in which, 
in all known examples, the scalar curvature of the information metric diverges at 
second order critical points, indicating a genuine (irremovable) singularity of the 
parameter manifold. This universal behavior can be related to the fact that near 
the critical points for CPTs, R is conjectured to be proportional to the correlation 
volume, which is not the case for QPTs, as we have stated. In fact, for CPTs, 
some components of the metric tensor are known to vanish at criticality (so that 
the geometry becomes ill defined, i.e, singular), but for QPTs, this is not the 
case. Exploring this direction further might be an useful exercise to undertake, 
particularly in the context of some known models in which classical analysis in 
the thermodynamic limit pinpoints the QPT as well [21], [2S]. Finally, in the 
absence of any singularity of the parameter manifold, it is interesting to ask how 
information geometry distinguishes between different phases of a theory that 
admits a QPT, in a coordinate independent manner. This issue will be addressed 
in a forthcoming publication. 
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